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ABSTRACT 


We consider the effect of the differential rotation of the disk of 
the Calaxy on magnetic field which penetrates the disk. The magnetic 
field will be progressively distorted from a potential (current-f rec.) form 
and will at some stage become unstable. We expect, from knowledge of 
solar flares, that an MHD instability, a resistive Instability, or u combin- 
ation of the two, will result in the release of the excess magnetic energy 
and that part of the released energy will be converted into heat. By 
estimating the energy release and the tate at which this process will 
occur and by assuming that this energy input is balanced by radiation, we 
obtain estimates of the parameters of the resulting plasma. It appears 
that this process alone can heat a galactic corona to temperatures of order 
10 6 K. 
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I. INTRODUCTION 


The diffuse soft x-ray (E < 1 keV) background has been attributed to 
a coronal-type plasma of T - 10^K which occupies a substantial fraction of 
the galactic volume (see Tanaka & Bleeker, 1977, for a review). Models 
for the creation of this high-temperature interstellar medium ( ISM) component 
through supernova explosions have been developed by Cox and Smith (1974), 
Smith (1977), and McKee and Ostriker (1977). Spitzer (1956), Shapiro and 
Field (1976), and Chevalier and Oegerle (1979) have analyzed the dynamics 
and energy balance of the hot ISM component; using the supernova heating 
me'hanism, they predict the formation of a corona extending several kilo- 
parsecs or more above the plan? of the Galaxy. We wish to suggest an 
alternative process for heating such a corona. 

Various lines of evidence indicate that the Galaxy is permeated by 
a magnetic field, possibly of complex structure. For the same reason that, 
if a magnetic field originates in the solar photosphere, it is energeti- 
cally favorable for the field to extend far out into the corona, we may 
infer that, if there is magnetic field in the disk of the galaxy, this 
field will also extend far outside the disk into what may be termed a 
"coronal" region. 

In the case of the sun, energy is released impulsively at coronal 
heights during flares (Svestka, 1976; Sturrock, 1979). It is generally 
agreed that the energy released during a flare is stored above the photo- 
sphere in the form of the free energy of a current-carrying magnetic field. 
The influence of differential photospheric motion on magnetic flux thread- 
ing the photosphere necessarily leads to non-potential field configurations 
and hence to free magnetic energy. 


2 



Given Che fact that the disk of the Galaxy is in a state of differ- 
ential rotation, and given the reasonable assumption that magnetic field 
threads the disk of the Galaxy and extends out into a galactic corona, 
the same Ingredients are present in the Galaxy as lead to flare activity 
in the case of the sun. The purpose of this article is to propose that 
galactic flares provide a substantial energy input to the Interstellar 
medium, and thus heat the coronal plasma responsible for the diffuse soft 
x-ray background. 
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II. ANALYSIS 


This analysis will be confined simply to order-of-magnltude estimates 
of the effects of flare activity on the galactic corona. The arguments 
and analysis will run parallel to an earlier similar discussion of solar 
flares (Sturrock, 1974). We consider what is Initially a closed magnetic- 
field structure of mean strength B (gauss) in the form of a flux tube. 
Because of the differential rotation of the Calaxy, both ends of such 
a flux tube are being twisted in the same sense, as viewed from the 
disk of the Galaxy, but in opposite senses with reference to the flux 
tube itself. If, at some time, a particular flux tube contains a potential 
field, then some Lime later the twisting of the field will lead to a 
helical configuration of the magnetic field. Unless the plasma stresses 
are large, the field will tend to adopt a force-free configuration. 

Theoretical studies indicate that a force-free field becomes I-fHD- 
unstable when the distortion from the current-free form reaches a critical 
level. Numerical studies (Barnes and Sturrock, 1972) indicate that, when 
the distortion exceeds the critical level, a closed magnetic-field con- 
figuration is likely to erupt to form an open configuration which neces- 
arily involves »:.e or more current sheets. These sheets can then reconnect 
by nonlinear processes (Petschek, 1964) arising from the tearing-mode 
instability (Furth et al., 1963), returning the field to its current-free 
form and, in the process, releasing a large fraction of the initial free 
energy associated with the coronal currents. 

Another somewhat simpler possibility is that a sheared magnetic-field 
configuration is at some stage subject to resistive instabilities which 
immediately return the field to the current-free form (Spicer, 1977). 
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Studies carried out during the recent Skylab Workshop on Solar Flares led 
participants to the view that both processes probably occur In the sun 
and are responsible for different types of solar flares (Sturrock, 1979). 

It is also possible (Sturrock and Uchida, 1979) that the heating of the 
normal corona may be due to the same processes occurring in magnetic fields 
of such configuration that the onset of Instability is non-explowfve 
rather than explosive (Sturrock, 1966). 

We suppose that, when a galactic flare occurs, energy in the amount 
1 2 

c jrj- B per unit volume is released in such a form as to heat the coronal 
plasma. The coefficient c will be the product of two terms, Ej and c^, 
where is the ratio of the magnetic free energy to the initial current- 
free magnetic energy, and e ^ is the fraction of the magnetic free energy 
which goes to heat the plasma. The calculations of Barnes and Sturrock 
(1972) indicate that flare action can occur after magnetic field lines 
are twisted by half a complete rotation and that, at that time, » 1. 

We make the simple assumption that, during an eruption, one half of the 
free energy of the force-free field survives as the free energy associated 
with the current sheet. We further assume that, when reconnection occurs 
in the current sheet, one half the energy goes into heat and the rest into 
mass motion. Hence we assume that e 0 *l/4 so that e*l/4. Assuming that 
the energy released during one flare is almost all radiated away by the 
time of the next flare, we see that, directly after a flare, 

3nkT m E -jr - B 2 (2.1) 

BIT 

where the gas is assumed to be almost pure hydrogen which is almost fully 

ionized and n is the density of either electrons or protons. 

If the plasma loses energy primarily by radiation, we may introduce 

the radiative-cooling time scale T defined by 

K 
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n • 3nkT x 1 

K 


( 2 . 2 ) 


-3 -1 

where n (erg cm s ) is the radiative rate. We adopt for n the form 

n - Tn 2 T Y (2.3) 

given by McKee and Cowie (1977) in the range 10^<T<10 2 ' 6 , where the 

— 18 2 

radiation rate is well represented by the above equation with T ■ 10 
and Y “ -3/5. 

The effect of differential rotation of the Galaxy is such that any 
element of a - ea exhibits an "orbital” rotation and, in addition, a "planet- 
ary" rotation which may be characterized by a vorticity given by 


<0 - 



(?. 4 ) 


where Q(r) is the angular velocity of the Galaxy at radius r. Hence a 

small flux tube which emerges from one location in the disk and returns 

to another location at about the same radius will be twisted at both ends 

in such a way that the flux tube is given a helical twist. We expect that 

the tube will become MHD unstable when the twist reaches a critical value. 

Numerical studies (Barnes and Sturrock, 1972) indicate that this value 

corresponds to half a complete rotation which would result from a rotation 

of n/2 at each end. Hence we may estimate the time between flares t (s) 

F 

from 


0) t_, ss tt/2 . (2.5) 

F 

We find from Mihalas and Routly (1968) that, over the range of radii 

3 kpc to 13 kpc, to has a value close to 30 km s *kpc * or 10 ^‘^s ^ . 

15 2 

Hence we obtain the estimate t_ = 10 

F 

If the corona plasma reaches an approximate steady state in which 
heating by flares is balanced by cooling by radiation, then 
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We obtain from (2.1), (2.2), and (2.3), the following estimates of the 
time-average /alues of the density and temperature: 


n . ID 9 3 B 16/13 t- 5/ ” 

T- 1C 4 ' 1 B 10/l3 T 5 ' 13 . 
which for the above estimate of T becomes 

„ - io 3 -V 6 ' 13 . 

T - ID 9 ' 9 B 10 ' 13 . 


(2.7) 


( 2 . 8 ) 


If the magnetic field, which contains the coronal plasma, extends out to 
a distance L (cm) from the disk, then the emission measure H, defined by 

H - n 2 L, (2.9) 


will have the value 


H - ID 7 ' 0 B 32/13 L . 


( 2 . 10 ) 



III. DISCUSSION 


Analysis of the soft x-ray background indicates that the plasma respon- 
sible for this emission probably has a complex thermal structure (see, 
e.g. Cash et al. f 1976; Bursteln et al . , 1977; Stem and Bovyer, 1979). 

Nevertheless, for the order-of-magnitude analysis presented here, we may 

5 • 8 

adopt as typical coronal parameters a temperature of ~ 10 * K and an emis- 
sion measure of * 10 16 ' 2 cm -5 ( 5 . 10~ 3 pc cm 6 ). We see from equation (2.8) 

that this temperature is consistent with the present model if the mean magnetic 

-5 . 3 

field strength in the local region of the galactic corona is ~ 10 ' gauss. 

For comparison, we note that Badhwar and Stephens (1977), on the basis of 
an analysis of the hydrostatic equilibrium of the gas-field system in the 
solar neighborhood, conclude that the magnetic field strength in the 
neighborhood of the sun is about (5-6) x 10 (> gauss, and that Ruzmajkin 
and Sokolov (1977), on the basis of an analysis of pulsar data, conclude 
that the local magnetic field strength is (2.1 ± 1.1) x 10 ^ gauss. Our 
theoretical estimate of the required magnetic field strength agrees with 
these two (observationally derived) estimates as well as they agree with 
each other. 

We now see also from equation (2.8) that, with the above estimate of 

the magnetic field strength, we expect the electron density in the galactic 

-3 0 -3 

corona to be about 10 " cm . We also see from equation (2.9) that the 

16 2 

inferred emission measure H b 10 ‘ is obtained if the corona extends out 

22 2 

to a distance L * 10 , i.e. 5kpc. 

It is necessary to compare the above estimate of L with the scale 
height h (cm) of plasma in the galactic corona. We ••-e 


kT 
mg * 


8 


(3.1) 


if the plasma is 


and note that the mean mass m (g) has the value 10 

mainly fully ionized hydrogen, and that the gravitational acceleration 
—2 -8 1 

g (cms ) is approximately 10 " in the range 1 kpc to 5 kpc of the 

galactic plane (Allen, 1973; p. 250), so that, numerically, 

h-10 16,3 T. (3.2) 

5.8 221 

We find that, for the estimated temperature T * 10 ’ , h * 10 . Hence, 

if the value of L is in fact set by the extent of the magnetic flux tubes 
which contain the plasma in the local region of the galactic corona, the 
scale height is sufficiently large to fill the flux tubes with plasma. 

If, on the other hand, the magnetic flux tubes extend well beyond 5 kpc, 
then the extent of the coronal plasma may be attributed to the scale height. 

If the latter is the correct interpretation, we may use (3.2) in (2.10) 
to obtain an expression for h in terms of B only. Then (as we have seen) 
our choice of B fits both T and H. However, our initial equations [(2.1) 
and (2.2)] are based on the assumption that n and T are uniform for ary 
flux tube, and this assumption must be relaxed if the scale height is less 
than the height of the flux tubes. 

The fact that the present model gives approximate order-or-magnitude 
agreement with observational data does not establish its validity; it simply 
indicates that more detailed investigation should be pursued. In particular, 
it is necessary to relax the "steady-state" assumption and to investigate 
the time development of heating, cooling and mass flow. The condition 
(2.1) is valid directly after a sudden release of energy, whereas condition 
(2.6) is valid for a long-term average of the coronal conditions. If 
were in fact short compared with x , then energy released during a flare 
would rapidly radiate away so that the time-average temperature would be 
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less Chan the Immediate post-flare temperature. If, on the other hand, 
t r were large compared to Tj. , then energy released during flares would accumu- 
late with the possible result that the plasma stress eventually exceeds 
the possible magnetic stress. This would result In an opening of the 
magnetic field with loss of the coronal plasma. 

If T (s) Is the time scale for the release of energy In a tiux tube, 

b 

chen (Sturrock, 1974) 


T 


E 



(3.3) 


where R (cm) is the radius of the flux tube under consideration and 

(cm s *) Is the reconnec 1. n rate. Following Petschek (1964), we esti- 
mate that 


v * 10 -1 v (3.4) 

R A 

where v is the Alfven speed. Since the mass density of fully ionize^ 

A 

hydrogen is 2nm, where m is the mean particle mass, 

v - B(8Trnm) _1/ 2 * 10 11,3 B n -1 ' 2 . (3.5) 

A 


Hence, using (2.7), we obtain 


t e - 10- 5 - 6 R If 5 ' 13 t - 5/26 


(3.6) 


or, for the value T - 10 


15.2 


10- 8 - 6 R B- 5/13 


(3.7) 


-5 3 

If we consider the value B = 10 , we find that, for R in the range 

10 23 ’ 3 to 10 22 ’ 2 (1 kpc to 5 kpc), T is in the range 10*^"^ to 10 33- ^. 


It appears, therefore, that T is comparable with Tj* This suggests that 
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Che corona wil . be approximately In a steady state so that the heating is 
more like the steady heating of an active region rather than the impulsive 
heating of a flare. 

We should also examine the ra"e of energy losr through thermal conduc- 
tion to the galactic disk, to compare with the energy loss rate by radia- 
tion. If T (s) is the time scale for cooling of the galactic plasma by 

thermal conduction, then T„ may be estimated from 

H 

3 nkTL- F T u , (3.8) 

H 

-2 -1 

where F (erg cm " s ) is the heat flux along the magnetic field. This 
may be estimated from 

F » tc T L'-IO" 6 T ?/2 L' 1 , 0.9) 

if we use the Spitzer (1962) form of the thermal conduction coefficient. 
Hence we obtain the estimate 


t - io- 9 - 4 l 2 „t- 5/2 . 

H 


(3.10) 


?? 2 -1 S ft 20 5 

For the values L ■ lO 4 "*” 4 ", n ■ 10 and T * 10 , we obtain T * 10 # 

H 

giving a cooling time of 10 l °y. This is much longer than the radiation 

cooling time, confirming that radiation is the dominant mechanism for energy 
loss. On the other hand, T„ is comparable with the age of the Galaxy, so 
that the energy flow from the corona to the galactic disk is sufficient to 
produce the "evaporation" (Neupert, 1968; Antiochos and Sturrock, 1978) 
necessar to provide the coronal plasma. 

It has been known for some time that there are high-velocity clouds 
(HVC) and intermediate-velocity clouds (IVC) of neutral hydrogen outside 
the galactic plane. Verschuer (1975), in reviewing this topic, suggests 
that there is more than one type of cloud and notes several possible inter- 
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pretatlons. Shapiro and Field (1976) have developed a "galactic fountain" 
model Involving an upward streaming hot gas which cools and condenses Intc 
a cool component. Our comparison between the galactic corona and the potu • 
flare plasma In the sun suggests that some of these clouds mav be the 
galactic analogues of "coronal rain" which often occurs In the loop promi- 
nence systems formed by some large solar flares (Tandberg-Hanssen, 1974). 

It has been shown by Goldsmith (1971) and by Antiochos (1976) that the 
latter may be understood as the result of thermal Instability In the hot 
dense plasma which forms in the coronal volume of an active region as the 
result of a flare. 

Following Field (1965) and Antiochos (1976), we note that the condi- 
tion for thermal instability Is that the wave number k (cm ) of the small- 
amplltui . »..vc should satisfy the inequality 


k 2 < 


k c -*"(? S - $)• 


(3.11) 


On using the form (2.3) for n and the Spltzer form for the thermal con- 
ductivity k , we find that 


. lrt 3 ,, *1/2 .1/2 _y/2 - 7/4 

k ■ 10 (2 - y) i n T 

c 

With the values of > and f appropriate to the range 10^ < T < 10^ * , 


becomes 


(3.12) 

this 


k 

c 


10 


-5.9 


n T - 41/20 

n T • 


(3.13) 


— 1 S 8 —1 20 8 

For the values n ■ 10 , f * 10 ' , we find that k ■ 10" * . Since 

c 

this estimate of k^* is considerably less than our estimate of L, we expect 
that thermal instability will lead to the formation of condensations 
with dimensions (along the magnetic field) of 200 pc or more. When formed. 
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these condensations should fall towards the galactic plane with speeds 

1/2 22 2 

approaching the free-fall speed (2 gL) . For L • 10 , this corres- 

7.2 -1 -1 

ponds to a velocity of 10 cm s or about 170 km s ; the Inferred 
speeds of high-velocity clouds are of this order of magnitude, and the 
Inferred speeds of intermediate-velocity clouds (which extend up to high 
galactic latitudes) are somewhat less. Hence it seems possible that some 
moving clouds are the galactic analogues of coronal rain, although 
Verschuur (1975) discredits the view that most clouds ave moving towards 
the galactic plane. 

Finally, we note that the proposed parameters for the plasma and 

magnetic field may be compared with data for the rotation of rtdio waves 

which enter the Calaxy from extragalactlc sources (Gardner and Whiteoak, 

1966). The rotation (radian) of the polarization of a -adio wave of 

wavelength X (meters) is normally written as 
m 

* “ Rm X 2 (3.14) 

where 

Rn" 10" 12 ' 6 fn B (| ds, (3.15) 

where B,| is the component of magnetic field along the propagation path. 
Hence we may estimate the maximum magnitude of from 

= 10 --0,6 n BL. (3.16) 

-5.3 - 3 22.2 

For the values B - 10 , n * 10 and L ■ 10 , this leads to a maximum 

value cf R^ of order 20. This estimate appears to be consistent with the 
upper limit of values of the rotation measure for high galactic latitudes 
(Gardner and Whiteoak, 1966, Fig. 9) . 
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